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Abstract 

We propose an efficient approach within the density-functional theory to determine the phonon 
structure of infinite solids under finite electric fields. We apply this approach to technological 
SrTiC>3, predicting many unusual properties — including a giant frequency shift, anomalous piezo- 
electric response, as well as striking dielectric tunability. A microscopic understanding for individ- 
ual phonon modes under finite electric fields is also provided. 

PACS numbers: 77.84.-s, 77.22.-d., 63.20.-e 
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It has been nearly 40 years since Worlock and Fleury discovered by means of Raman 
scattering that under finite electric fields the transverse-optical (TO) phonon in SrTiC>3,Q] 
as well as in KTa03,|2| exhibits a striking shift in frequency — by more than 400% of its 
zero-field value — when a very small field of 12 kV/cm is applied. Q] Despite its obvious 
importance, this giant field-induced shift of phonon frequency has not thus far been inde- 
pendently investigated and/or understood via the powerful density-functional theory (DFT) 
which, on the other hand, has been amply applied in the past three decades on a variety of 
other material properties. The reason lies in the difficulty in handling finite electric fields 
within the direct DFT approach, a topic that has recently drawn a great deal of interest 

and progress. □□□□□□□ 

SrTiC>3 is an incipient ferroelectric that is approaching, but does not in fact undergo, a 
paraelectric-to-ferroelectric phase transformation when lowering temperature, and is often 
considered as an "inert" perovskite. jloT| However, being on the edge of phonon softening, 
SrTiC>3 shows unusual sensitivity on external electric fields that are able to drive easily 
the occurrence of the phase transition, and thus likely lead to an anomaly in the field- 
induced structural response. Given their technological importance, this and other similar 
ferroelectric incipients (in particular, when they are under finite electric fields) remain, 
nevertheless, poorly understood. For instance, it is not clear how different phonon modes 
may respond differently to the electric fields in terms of their frequencies and in terms of 
their atomic participations (that is, the phonon properties). Furthermore, little is known 
on how atomic displacements driven by external fields will couple with strain in incipients 
(i.e., the mechanical properties), and whether SrTiC>3 in finite electric fields can be a good 
piezoelectric, in addition to being a good dielectric. Finally, a direct DFT modelling of the 
"dielectric tunability" by finite electric fields in incipients (i.e., the dielectric properties) is 
still at its early stage and thus much needed. 

The purpose of this Letter is threefold: (1) to show that it is possible to devise a direct 
density-functional approach to determine and understand the phonon spectrum of infinite 
solids under finite electric fields; (2) to demonstrate that the approach works robust for 
incipient SrTiC>3, and in particular, the theoretical prediction of field- induced frequency 
shift is in excellent agreement with the result observed in experiments; (3) to present other 
unusual properties in SrTi03 that occur "only" under finite electric fields, that is, a gigantic 
piezoelectric response as well as a superior field-induced dielectric tunability. 
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To determine phonon structure, namely vibration frequencies and phonon eigenmodes, of 
a solid under finite electric fields, we start with free energy F = U — E- P, where E (assumed 
here to be along the pseudo-cubic [001] direction) and P are the screened macroscopic field 
and polarization per cell, respectively. This yields force-constant matrix elements as 

d 2 F d 2 U dZjJi) 

, I 1 ) 



dr ia dr jf3 dr ia dr jf3 dr jf3 

where r ia is the position of atom i along the a cartesian direction, and Z*(i) is the Born 
effective charge tensor of this atom. It should be emphasized that each term in the right- 
hand side of Eq.((TJ) must be evaluated at field-induced new atomic-positions {R new } and 
cell-shape strain {??new}, instead of the zero-field {Ro} and {f]o}. 

Determination of the force-constant matrix of Eq.Q is done in two steps. First, {Rn™ 1 
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and {^new} are obtained using a newly-developed method that was described in Ref. 
Second, the force- const ant matrix is evaluated as follow. The first term in Eq.(fTj). which is 
associated with the zero-field internal energy U but at displaced atomic positions, is ready 



to be computed using standard density-functional perturbation theory (DFPT). 13j, [14] The 
second term requires, however, a perturbation theory of the macroscopic polarization P with 
respect to the atomic displacements, and is rather cumbersome to deal within the DFPT. 
Here, we numerically find (by performing finite-difference calculations) that, even for an 
electric field as large as 0.5xl0 8 V/m, the second term in Eq.(jT|) turns out to be less than 
0.5% of the first term for all diagonal elements, and thus contributes practically very little 
to the dynamic matrix. This particular finding allows us to neglect the insignificant second 
term. The resulting simplicity and efficiency in the present approach is precisely what we 
pursue, and allows the determination of the phonon structure under finite electric fields to 
be routinely performed within the DFT. This approach works in general for phonon at any 
wave-vector q, whereas in the following we focus on the phonon modes at the zone center 
since they can be directly compared with Raman experiments. 

0,0 

It should be pointed out that, apart from the (externally-applied) field E, a long- 
wavelength longitudinal-optical (LO) phonon in polar materials will generate a second 
macroscopic field Ei with its direction depending on the phonon wave-vector q. Here we 
choose to study a configuration that likely leads to most pronounced effects that are to be 
caused by the applied field E, namely, when E and q are perpendicular, with q chosen along 
the pseudo-cubic [010] direction. 
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Obviously, quantum fluctuation effects in SrTiOs at low temperature, which are respon- 
sible for its incipient behavior, can not be described by zero-temperature DFT calculations 
and are not the focus of this work. On the other hand, the DFT is able to predict quan- 
titatively the electric- field effects in SrTiOs, since these effects are mainly determined by 
the zero-field frequency (i.e., the energy surface) of the ongoing soft mode (assuming that 
this frequency is correctly reproduced by DFT). Furthermore, anti- ferroelectric rotation in 
low-temperature SrTiOs is not considered here, as it affects mainly the zone-boundary (not 
the zone-center) phonons. 

n 

Technically, our calculations are carried out using the ABINIT code.pjfl modified to 
perform finite-field geometry optimization. A sufficient energy cutoff of 110 Ry is used 
to ensure a satisfactory convergence. The zero-field lowest TO phonon, that turns soft 
in normal ferroelectrics, is found to be still stable and at the frequency of 13.2 cm -1 for 
cubic SrTiOs with the lattice constant 3.892 A. This frequency is in good agreement with 
experimental value of 11 cm -1 [J, and in fact, an agreement is found for all modes across 
the phonon spectrum. 

FigO] shows the phonon spectrum under an electric field of 0.22xl0 8 V/m, compared to 
the zero- field case. The most striking results are: (1) Upon the electric field, the degeneracy 
is completely lifted for all phonon modes, and in particular, for the zero-field "silent" modes 
(being triply-degenerate with r 25 symmetry in cubic SrTiOs). Being non-polar with only 
oxygen amplitude, these silent modes were commonly thought to split into one singlet and 
one doublet by electric fields jlfil ]. (2) The lowest two modes in the lower panel of FigJH 
which split from a doublet TO phonon in the upper panel, undergo giant shifts in frequency 
from 13.2 cm -1 at zero field to 48.7 and 92.2 cm -1 at the considered magnitude of field. In 
other words, the frequencies of these modes are altered by the applied electric field by more 
than 400% and 800% (!), respectively. (3) Another mode at zero-field frequency of 163.2 
cm -1 (notice that this is not the ongoing soft mode) displays also a remarkable frequency 
shift by as much as 26.6 cm -1 . (4) On the other hand, modes at 447.0 and 772.4 cm -1 are 
found to be least affected (if any). 

In order to understand the microscopic origin of why some modes, not others, are drasti- 
cally altered by electric fields, and furthermore, to provide valuable guide for understanding 
future experimental results, Table |l] gives the vibration eigenvectors of the phonon modes. 
Surprisingly, the second mode in this Table — namely, the mode displaying the largest field- 
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induced frequency shift, vibrating in parallel with the E direction, and its frequency to be 
referred to as uu — turns out to have virtually no amplitude from Ti, instead it is predom- 
inantly from strontium and oxygen atoms. This is in contrast with the zero-field situation 
where Ti contribution is large (~ 0.31) for the corresponding mode. On the other hand, 
another (also z-axis polarized) mode at 189.8 cm -1 shows an opposite trait, with a heavy Ti 
involvement. Also interestingly, it can be seen in Table |l] that a large field-induced frequency 
shift occurs not only to those modes that vibrate parallel to E, as it should (for instance, 
see the modes at 92.2, 189.8, and 549.9 cm -1 ) — and in fact, the applied field also produces 
a remarkable shift to the mode that is transverse to the field direction (see, the first mode 
in Table HJ and the frequency of this mode will be denoted as u± hereafter). Meanwhile, the 
modes that exhibit the least frequency shift (e.g., modes at 447.0 cm -1 and at 772.4 cm -1 ) 
are found to associate with vibrations along the [010] direction. This can be explained by 
the fact that these modes originate from the E x field due to the LO/TO splitting and thus 
subject to little effects from the applied E field. 

We now wonder how the phonon frequencies may vary with different field-strengths, 
and the results are depicted in FigJ2k for the u)u and u>± modes. Intriguingly, our results 
show that the field- induced increase in frequency (i.e., so-called "phonon hardening") is 
characterized by two drastically different regimes, separated by a critical field E c and equal 
to ~0.5xl0 8 V/m. Below E c , un and uj± increase dramatically and strong nonlinearly — 
with their difference u\\ — u>± first enlarged and then saturated at E c — as the field strength 
increases. On the other hand and above E c , u\\ and u± show almost a linear field-dependence, 
and interestingly, the frequency-vs-field slopes are nearly equal for these two modes, whereas 
they vibrate along completely different directions. Fitting the two curves in Fig|2^ at high 
field-strength (E > 0.25 x 10 8 V/m) to an analytical expression uj 2 (E) + uo 2 (E = 0) = oE 1 
(derived from the Ginzburg-Landau-Devonshire theory [lj) yields, for the u>\\ and u± modes, 
an exponent 7 of 0.31 and 0.56, respectively. Our direct DFT results thus demonstrate 
that the exponent 7 not only significantly varies from mode to mode, but also considerably 
deviates from the value of 0.66 from the phenomenological theory. |l7j] 

The peculiar behavior on the field-dependence of the mode frequency in Fig|2t is found 
to arise from the atomic displacements and cell-shape changes induced by the electric fields, 
shown in Fig|2b- It can be clearly seen that below the critical E c field the dramatic and 
nonlinear displacements of O atoms (with respect to Ti), as well as the sharp increase of the 
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c/a strain, correlate well with the phonon responses in FigJ2K Notably, two non-equivalent 
oxygen atoms (namely, Oi and O3) respond remarkably similar to the applied electric field, 
whereas the (relative) displacement of Sr is almost negligible. These results thus show 
(not surprisingly) that the field effect is primarily to move Ti with respect to the rather 
rigid oxygen octahedral, while Sr contributes by having the "right" atomic size to form an 
incipient (Indeed, the low-frequency modes in Table |1J which matter most for the incipient 
behavior, all show a significant Sr participation). 

Our theoretical attempt on phonon under finite electric fields yields results that quan- 
titatively agree with experiments. For example, besides a sharp Raman peak that results 
from the cou phonon, one additional peak, despite its weak intensity, was clearly observed at 
a lower frequency. |l| This additional peak is consistent with the existence of the ujj_ mode as 
predicted in our calculations. Furthermore, our theoretical u\\ frequency at the 0.012xl0 8 
V/m field is determined to be 42.7 cm -1 , in excellent agreement with the experimental value 
of 45 cm -1 , measured under the same field strength and at temperature 8K.Q] 

Finally, we would like to examine other technological properties — namely, piezoelectric 
633 — e3i coefficient (which measures the piezoelectric response to a tetragonal strain when 
volume is preserved) and static dielectric 5 33 and en constants — of SrTi0 3 that is placed 
under finite electric fields. The 633 — e^\ coefficient is null at zero field (see, FigJ2K)- This 
coefficient becomes, however, strikingly large — being as much as 19 C/m 2 under a small field 
of 0.012x10 s V/m. Note that this magnitude of piezoelectric response in SrTiC>3 (which is 
predicted here for the first time, to our knowledge) exceeds far more than the value of 

n 

4.08 C/m 2 in PbTi0 3 [18], and even than the intrinsic theoretical 633 value of 9.1 C/m 2 
in PMN-PT ^j] . We attribute this exceptional response to the field-induced occurrence 
of the paraelectric-to- ferroelectric phase transition in SrTiOs, which does not happen in 
PbTiOs and in PMN-PT. Interestingly, the 633 — e3i coefficient in Fig|2K first increases, 
then decreases, and eventually levels off at a nonzero value of ~ 6 C/m 2 . Regarding the 
dielectric properties, the £33 dielectric constant is determined to be 12215, 1120, 188, 117, 
77, and 61 for the six considered field strengths at 0, 0.012, 0.22, 0.50, 1.0, and 1.5 xlO 8 
V/m, respectively. This tells us that the £33 constant decreases dramatically by more than 
one order of magnitude when going from zero field to a very small field of 0. 012 x10 s V/m, 
demonstrating the existence of an incredible field-induced dielectric tunability! This giant 
tunability is found also true for the £%i component, being 12215, 7123, 866, 485, 337, 264 at 
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the field strengths given above. Last (but not least), we want to point out that our direct 
DFT calculations of phonon frequencies and (independently) dielectric constants prove that 
the Lyddane-Sachs- Teller (LST) relationship [20] is fully valid under finite electric fields. In 
fact, the static dielectric constants obtained from direct DFPT calculations and those from 
the LST relationship (using phonon frequencies) are found to differ by less than 0.1%. 

In summary, a practical and efficient DFT approach that allows to determine the phonon 
structure of solids in finite electric fields was developed, and applied, to study the electric- 
field effects on the zone-center modes of incipient SrTiC>3, yielding a good agreement with ex- 
periment. Our results confirm quantitatively the striking field- induced frequency shifts that 
have been experimentally observed — but not explained by direct first-principles theory — for 
nearly four decades. Our work further provides a microscopic understanding of individual 
phonon mode under finite fields. We also find: (i) A giant phonon hardening (as the fre- 
quency increases by several fold with respect to the zero-field value) takes place, not merely 
for the uj\\ mode that vibrates in parallel with — but also for the u±_ mode that, on the 
other hand, vibrates perpendicularly to — the direction of the applied electric field, (ii) After 
hardening, the mode showing the largest frequency shift (i.e., the uj\\ mode) has a surpris- 
ingly small Ti-atom participation, (iii) The mode at the zero-field frequency of 163.2 cm -1 , 
whereas not being the ongoing soft mode, displays also a remarkable response to the im- 
posed electric field. Modes at 447 and 772 cm -1 are found to be least affected by the E field, 
(iv) The field dependence of the soft-mode frequency exhibits a crossover at E c =0.5xl0 8 
V/m — below which u\\ and u± upshift sharply and nonlinearly, but differently. Above E c , 
the responses of these two modes are interestingly similar with a nearly equal frequency-vs- 
field slope. This dependence is found to result from the field effects on atomic movements 
and on the tetragonal c/a strain, (v) A small electric field of 12 kV/cm is able to decrease 
the £33 constant by one order of magnitude, showing an extraordinary dielectric tunability. 
(vi) An unusual piezoelectric response is predicted to exist in SrTi0 3 when placed under 
finite fields. We expect that the present approach will have widespread applications for 
determining phonon structure and dielectric tunability in solids under finite electric fields. 
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FIG. 1: Mode frequency and degeneracy of SrTi03 under zero field (the upper panel) and under 
a field of 0.22xl0 8 V/m (the lower panel). Zero-frequency acoustic modes are not shown. Dotted 
lines indicate the evolution of selected phonon modes. 

FIG. 2: Dependence of (a) the oj\\ and uj± frequencies, and (b) the z-direction displacements of Sr, 
Oi and O3 atoms (using the left vertical-axis) and the c/a strain (using the right vertical-axis), as 
a function of field strength. Inset in part (a) is the piezoelectric 633 — e^i coefficient. In part (b), 
Ti is chosen as the reference for measuring the z-axis displacements of other atoms. Symbols are 
calculation results; lines are guide for eyes. 

TABLE I: Frequency u and mode eigenvector £ T of the zone-center phonons in SrTi03 under a 
field of 0.22xl0 8 V/m, as well as frequency shift Aw = 00(E) — lo(E = 0). One eigenvector has 
three vibration components £ x , £ y and £ 2 (as specified by index r); each component is described 
using the notation £ T = (us T ,UTi,uo 1 ,uo 2 ,uo 3 ) with five amplitudes from Sr, Ti, Oq, O2, and O3 
atoms, respectively. O3 is located beneath the Ti atom. Those components not given in the table 
are zero. The x, y, and z-ax.es are, respectively, pointing along the pseudo-cubic [100], [010], and 
[001] directions. 



uj (cm v 


Aw (cm 1 


eigenvectors 


48.7 


35.5 




{ 0.43, 0.27, -0.43, -0.53, -0.52) 


92.2 


79.0 




;-0.57, -0.04, 0.51, 0.51, 0.40) 


149.9 


1.2 


Zy= 


;-0.71, 0.60, 0.23, 0.16, 0.24) 


166.3 


3.1 


6r = 


(-0.58, 0.82, -0.02, -0.03, 0.01) 


189.8 


26.6 




;-0.44, 0.86, -0.17, -0.17, -0.11) 


226.0 


2.2 


Hy= 


;-o.oi, 0.01, 0.71, 0.00, -0.71) 


226.4 


2.6 


6r= 


(-0.01, 0.02, 0.00, 0.70, -0.71) 


227.0 


3.2 




; 0.00, 0.00, -0.71, 0.71, 0.00) 


447.0 


0.0 


iy= 


(-0.05, -0.43, 0.59, -0.34, 0.59) 


542.3 


1.2 




( 0.04, 0.02, -0.85, 0.37, 0.37) 


549.9 


8.8 




; 0.03, 0.04, 0.36, 0.36, -0.86) 


772.4 


0.6 




(-0.12, -0.45, 0.08, 0.88, 0.08) 
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